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NOTATION 

Source  amplitude 

Exponent  of  density  stratification 
Phase  velocity 
Acceleration  due  to  gravity 
Magnitude  of  group  velocity 
Components  of  group  velocity 

Hankel  function  of  the  first  kind  of  order  v 
Magnitude  of  wave  number 

Modified  Bessel  function  of  the  second  kind 
of  order  v 
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CD 

Valsala  frequency 
p2  =  _  M^ 

Total,  perturbation, and  mean  pressure,  respectively 

Source  distribution 

Time 

Components  of  perturbation  velocity 
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Components  of  total  velocity 
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INTRODUCTION 

Love  (see  Lamb  19^5 was  the  first  to  derive  the  linearized 
equations  governing  small  oscillations  about  the  state  of  hydro¬ 
static  equilibrium  In  a  continuously  stratified  fluid.  The  propa¬ 
gation  of  free  elementary  waves  In  continuously  stratified  fluids 
were  studied  by  Ylh  (i960),  Eckart  (i960),  Yanowltch  (1962)  and 
Tolstoy  (1963)-  Some  of  these  authors  also  considered  the  effects 
of  compressibility  and  rotation. 

In  this  report  the  Internal  waves  generated  by  a  time  de¬ 
pendent  local  disturbance  In  an  Incompressible  continuously  strat¬ 
ified  fluid  are  examined.  Following  the  Love  theory  of  wave  motion 
In  heterogeneous  liquids,  the  motion  Is  regarded  as  a  small  per¬ 
turbation  about  the  state  of  hydrostatic  equilibrium.  For  sim¬ 
plicity,  the  equilibrium  density  profile  Is  taken  to  be  expo¬ 
nential  whereby  the  coefficients  of  the  governing  equations  be¬ 
come  constants.  The  complications  of  multiple  reflections  from 
physical  boundaries  are  avoided  by  considering  a  fluid  of  Infinite 
extent.  The  disturbance  Is  taken  to  be  two-dimensional  and  time 
dependent  but  non-moving. 

THE  GOVERNING  EQUATIONS 

Consider  the  equations  governing  the  two-dimensional  motion 
of  an  Incompressible,  continuously  stratified  fluid  under  gravity. 
Let  (x,y)  be  a  Cartesian  coordinate  system  oriented  with  gravity 
pointing  In  the  negative  y  direction,  and  (u,v)  be  the  velocity 
components  In  the  x  and  y  directions  respectively. 
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Slnce  the  fluid  le  incompressible,  then 


du  ^ 
bx  by 


27rq(x,  y,  t) 


[1] 


where  27rq  denotes  the  volume  rate  of  outflow  from  a  source  dis¬ 
tribution  in  the  fluid.  Conservation  of  mass  demands  that 


dt 


+ 


P 


Su  ^ ' 
,Sx  by  I 


hp 


^  = 


27rq(x,y,  t)p^(x,y,t) 


[2] 


where  p  is  the  density  field  and  p^  is  the  density  of  the  fluid 
being  given  out  (in  this  case,  p^  is  left  open)  or  taken  in  (in 
this  case,  p^  must  be  identically  equal  to  the  ambient  density) 
by  the  source  distribution.  Subtracting  [l]  from  [2]  yields  the 
following  'equation  of  state'; 


St 


+  u 


Sx 


+  V 


^y 


2Trq(p^-  p)H(q) 


[3] 


where  H  denotes  the  Heaviside  function.  If  p^  is  controlled 
such  that  it  is  always  equal  to  the  ambient  density,  then  Equa¬ 
tion  [3]  simplifies  to 
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ie. 

St 


+  u 


Sx 


+  V 


^  = 


0 


[^] 


In  v;hat  follows,  this  is  assumed  to  be  always  the  case. 

Equations  [l]  and  [4],  together  with  the  Euler  equations 


Su  Su  „  Su  1  Sp 

St  Sx  Sy  ~  ■  p  Sx 


[5] 


Sv  ,  Sv  ,  ,  Sv  1  Sp 

St  ^  Sx  Sy  p  Sy 


g 


[6] 


where  p  denotes  the  pressure  and  g,  the  acceleration  due  to 
gravity,  form  a  system  of  four  equations  for  the  four  unknowns 
(u,v,p,p). 

If  the  motion  Is  weak.  It  may  be  regarded  as  a  small 
perturbation  about  the  state  of  hydrostatic  equilibrium  and, 
as  a  first  approximation,  second  order  effects  may  be  neglected. 
Denoting  the  perturbations  by  (u',  v',  p',  p’ ),  the  hydrostatic 
pressure  by  p  and  the  equilibrium  density  by  p,  Equations  [l], 
[4],  [5]  and  [6],  on  neglecting  second  order  terms,  become 


4 
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du ’  dv  * 
dx  d  y 


2v’q 


[7] 


0 


-  du'  dp ' 

P  ^  ■  dT 


[9] 


p 


dv'  _  dp ' 
dt  ~  ~  dy 


gp' 


[10] 


The  coefficients  are  seen  to  be  functions  of  y  only  and  the 
resulting  system  is  linear.  From  this  system  of  equations,  the 
following  equation  for  v'  may  be  uncoupled: 


dy' 


(pV  ) 


dx' 


(pV  )  + 


Ml 

g 


_d_ 

dy 


.  2N  dN  -  ,v 

(pv)  +  -g- d7  P''  V 


+  N' 


dx' 


(pv'  )  =  2v- 


p^ 

dt^ 


[11] 


where  N  is  the  Valsala  frequency  defined  by 
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_r,  _ 

This  Important  parameter  serves  to  characterize  the  capability 
of  the  medium  as  an  internal  wave  carrier. 

In  what  follows,  attention  is  focused  entirely  on  the 
vertical  velocity  component  only,.  It  is  convenient  to  transform 
it  to  a  new  variable  V  according  to 


V  =  VT  V 


[13] 


The  Equation  for  V  is 


0  +  ^  + 


N  ^ 
g  by 


4g2 


V 


+  N‘ 


_d_ 


\ 


dt^ 


I 


[14] 


For  an  exponential  density  profile,  viz. 


P 


e-2by 


where  and  b  are  constants.  Equation  [l4]  simplifies  to 
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o^V  ^ 


b^V  )  +  N' 


2  dfV 

dx^ 


[16] 


and  the  Valsala  frequency  N  =  V  2gb  Is  constant  throughout  the 
fluid  which  means  tliat  the  fluid  Is  homogeneous  In  character  as 
a  wave  carrier  but  anisotropic . 

ELEMENTARY  WAVES 

In  this  section, certain  results  concerning  the  behaviour 
of  elementary  waves  In  a  stratified  fluid  are  discussed .  For  a 
more  general  treatment  In  this  respect,  ti‘;e  reader  Is  referred 
to  Eckart  (i960). 

The  homogeneous  character  of  the  fluid  suggests  that  plane 
waves  of  the  form 


1  (ax  +  jSy  -  cut ) 

V  =  A  e 


where  A  is  a  complex  constant,  a,  and  co  are  real  constants, 
are  possible  homogeneous  solutions  of  Equation  [I6].  Substitution 
into  Equation  [16]  reveals  that  this  is  Indeed  so  provided  that 
the  following  characteristic  equation  Is  satisfied,  viz 
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CJD^  (a^  +  =  0  [1?] 

In  the  first  place,  a  necessary  condition  for  Equation  [l?] 
to  have  real  solutions  Is  that  oi  must  be  less  than  the  Valsala 
frequency,  l.e. 


2 

—  -  1  =  >  0 

a 

CD 


Since  the  Valsala  frequency  measures  the  stability  or  'stiffness' 
of  the  fluid,  the  above  condition  Implies  that  higher  frequency 
waves  require  a  stiffer  medium^  a  result  which  Is  In  accord  with 
•Intuition.  For  each  q  that  is  less  than  the  Valsala  frequency. 
Equation  [17]  Is  satisfied  by  all  points  lying  on  an  hyperbola 
of  two  sheets  in  the  wave  number  plane.  In  Figure  1  the  hyper¬ 
bolas  in  the  first  quadrant  for  various  cd  are  shown.  The  family 
of  curves  Is  symmetric  about  both  the  vertical  and  horizontal 
axes.  It  Is  seen  that,  for  a  given  co/N,  the  magnitude,  k,  of 
the  wave  number  depends  on  Its  Inclination  from  the  horizontal 
axis,  0,  which  cannot  exceed  tan"^  M. 

The  phase  velocity  or.  In  other  words,  the  velocity  of 
propagation  of  the  Individual  crests  and  troughs  has  the  same 
direction  as  the  wave  number  and  Its  magnitude,  c,  is  given  by 
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c 


N 

b 


cos^S 


[18] 


Figure  2  is  a  polar  plot  of  c  against  9  for  various  values 
of  00.  The  fact  that  c  is  a  function  of  0  demonstrates  the 
anisotropic  character  of  the  medium. 

The  group  velocity,  according  to  which  energy  is  propagated, 
is  given  by 


(ax.ay)  =  [|S 


8/3  ; 


[19] 


The  direction  of  the  group  velocity,  then,  is  normal  to  that 
hyperbola  in  which  the  wave  number  is  found  at  (a,/3)  and  is,  in 
general,  different  from  the  direction  of  tlie  phase  velocity. 

This  behaviour  is  characteristic  of  wave  motions  in  anisotropic 
media  (see  Llghthlll  i960).  At  6  =  0  the  two  velocities  have 
the  same  direction,  as  6  Increases  the  difference  in  direction 
increases  and  tends  to  Tr/2  as  0  tends  to  tan“^  M.  The  magnitude 
of  the  group  velocity,  G,  is  found  to  be 
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Flgure  3  Is  a  polar  plot  of  G  against  Q  for  various  values  of  cu. 

It  Is  noted,  that  the  amplitude  of  the  actual  velocity,  v'. 
Is  not  constant  but  varies  with  y  as  l/V^  .  It  follows  that 
the  first  order  kinetic  energy  Is  conserved. 

INTERNAL  WAVES  GENERATED  BY  AN  OSCILLATING  SINGULARITY 

Consider  the  waves  generated  by  a  source  of  constant  ampli¬ 
tude  A  and  simple  harmon]  c  time  dependency  of  frequency  oo  located 
at  the  origin.  After  a  sufficient  time  has  elapsed  for  the 
transients  to  die  away,  the  dependent  variables  would  also  ac¬ 
quire  the  same  tim.e  dependency.  Therefore, 


-loot 

V(x,y,t)  =  ct(x,y)e 


[21] 


The  equation  governing  the  amplitude  function  cl(x,y)  Is 


^  ^  -  b^iJ>  =  GttA  V~^  5(x)[6(y)  -  b6(y)]  [22] 

dy^  dx^  ° 


where  6  denotes  the  Dirac  delta  function  and  a  dot  denotes  a 
differentiation . 

First,  consider  the  Green's  function,  <]),  of  Equation  [22]. 
By  definition,  4>  satisfies  the  following  equation; 


HYDRONAUTICS,  Incorporated 


-10- 


! 


^-1  _  _  ^20  =  6(x)6('y) 

dy®  dx® 


The  two  dimensional  Fourier  Integral  representation  of  <t>  Is  easily 
found  to  be 


-hoo  -hoo 


i(ax+  jSy) 


(t>(x,y)  = 


(2Tr)= 


b^) 


dadjS 


[24] 


Consider  first  the  case  where  cs  is  larger  than  the  Vaisala 
frequency,  thus  Is  real  and  negative.  Let  =  -  P®  where  P 
Is  real  and  positive.  Rewrite  Equation  [24]  as 


4) 


l(ax+  jSy) 
e 


a®P®  +  (/3®+  b®) 


-00  -00 


dad^ 


[25] 


Integrate  first  with  respect  to  a.  The  Integrand  has 
simple  poles  at 


d 


+ 


P 


[26] 
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The  contour  is  closed  by  a  large  semi-circle  In  the  upper  or 
lower  half-plane  depending  on  whether  x  Is  positive  or  negative 
and  the  Integral  evaluated  by  the  method  of  residues;  the 
contribution  from  the  semi-circle  vanishes  as  Its  radius  tends 
to  Infinity.  From  here  on,  attention  Is  focussed  on  x  >  0  only. 
For  X  >  0 


But 


-00 


-  f  V  b= 

Vi3^+  b^ 


cos  /3y  dj3 


cos  j3y  d/3  =  K 


b|^+  y^ 

P 


(see  Watson  1958)>  where  denotes  the  modified  Bessel  func¬ 
tion  of  the  second  kind  of  order  v.  Therefore, 
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<t> 


[27] 


Next,  consider  the  case  where  u)  is  less  than  the  Vaisala 
frequency,  thus  Is  real  and  positive.  The  poles  now  lie  on 
the  real  a  axis  and  the  integral  Is  indeterminate.  Under  such 
circumstances  it  is  usual  to  deform  the  path  of  integration 
slightly  such  that  it  goes  around  Instead  of  through  the  poles. 
How  the  path  should  be  deformed  depends  on  additional  physical 
information.  In  the  present  problem  it  is  physically  obvious 
that  must  be  symmetric  in  x  and  that  all  internal  waves  are 
outgoing  waves  from  the  source  at  the  origin.  With  this  con¬ 
siderations  in  mind  it  is  found  that  the  path  of  integration 
should  be  displaced  above  the  negative  pole  and  below  the  posi¬ 
tive  pole  by  small  semi-circles.  The  Integral  is  then  evaluated 
as  in  the  previous  case  to  give 


<t> 


+00 


*00 


i(^y  +  I 


djS 


1  I  V  ^^+b^ 


1 

2tM 


o 


Vi3^+  b® 


cos  jSy  d^ 
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But 


1 


cos  /3y  d/3  =  ^ 


U/? 

iw^ 

\  yC 


r 


(see  Watson  1958),  where  denotes  the  Hankel  function  of  the 

first  kind  of  order  v.  Therefore, 


4m 


H 


(1) 


[28] 


The  motion  generated  by  the  oscillating  source  can  be  ex¬ 
pressed  In  terms  of  <t>  as 


V  (x,y,t) 


-lODt 

e 


'2TrA 


I  8/^ 

8y 


[291 


First  examine  the  case  where  cd  Is  greater  than  the  Valsala 
frequency.  Substituting  Equation  [27]  Into  [29]  the  following 
equation  for  v'  Is  obtained 
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V' 


+  K 

o 


' 

b 

I 


[30] 


As  the  stratification  tends  to  zero,  l.e.  as  b  -►  0 


V' 


-la^t 

e 


A 


x^+ 


[31] 


This  limit  Is  the  motion  generated  by  an  oscillating  source 

In  a  homogeneous  fluid.  For  finite  b,  the  standing  wave  character 

of  the  motion  Is  preserved  because  both  K  and  Ki  are  real  for 

o 

real  arguments.  However,  the  purely  antisymmetric  character  of 
the  vertical  velocity  component  with  respect  to  y  Is  destroyed. 

It  Is  noted  that  no  progressive  internal  wave  Is  generated.  This 
is  in  accord  with  the  result  of  the  last  section  which  excludes 
the  possibility  of  progressive  waves  whose  frequency  Is  greater 
than  the  Valsala  frequency. 

Next  consider  the  case  where  cd  is  less  than  the  Valsala 
frequency.  Substitution  of  Equation  [28]  Into  [29]  results  in 
the  following  equation  for  v': 


‘3 


a 
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V' 


-lout 

e 


by 

TrAb  e 
2M 


[32] 


The  solution  behaves  vei-y  differently  depending  on  whether 
x^/M^  -  y^  is  positive  or  negative.  Consider  first  the  region 
which  satisfies  the  Inequality 


y^  <  0 


Rewrite  Equation  [32]  as 
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Uslng  the  following  relation,  (see  Watson  1958) 


ivTTl  ,  . 

K^(0  =  l^rl  e  (iC) 


where  ^  is  real  and  positive.  Equation  [33]  becomes 


It  Is  seen  from  the  behaviour  of  K  and  Kj  that  the  motion  in 

o 

this  region  is  a  standing  wave  also,  but,  here,  the  amplitude  of 
the  motion  is  singular  on  the  lines  y®  =  x^/M^  Instead  of  Just 
at  the  origin  where  the  source  is  located. 

In  the  region 


-  >  0 
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the  arguments  of  the  Hankel  functions  appearing  In  Equation  [32] 
are  real .  From  the  known  behaviours  of  the  Hankel  function  for 
real  arguments  It  Is  seen  that  internal  gravity  waves  are  present 
In  this  region. 

This  peculiar  pattern  of  motion  is  related  to  the  behaviour 
of  elementary  waves  discussed  in  the  previous  section.  The  period! 
source  may  be  regarded  as  an  Integral  of  elementary  waves  of  all 
wave  numbers  and  of  a  given  frequency  oi.  It  was  seen  that  only 
those  waves  whose  number  lie  on  a  particular  hyperbola  are  pro¬ 
gressive.  The  asymptotes  of  the  hyperbola  make  an  angle  of 
tan~^  M  with  the  horizontal  so  that  the  direction  of  the  wave 
number  cannot  exceed  tan"^  M.  However,  the  group  velocity  or 
the  velocity  of  propagation  Is  normal  to  the  hyperbola  so  that 
the  region  In  whj oh  progressive  waves  are  found  are  sectors 
bounded  by  the  lines  y®  =  x^/M^  as  shown  In  Figure  4.  Except  for 
the  small  wave  number  range  the  group  velocities  of  the  elementary 
waves  are  approximately  In  the  direction  of  either  one  of  these 
lines.  Thus,  there  Is  a  large  concentration  of  energy  in  them. 

INTERNAL  WAVES  GENERATED  BY  AN  INSTANTANEOUS  SOURCE 

In  this  section,  the  transient  problem  of  the  motion  due  to 
a  source  at  the  origin  with  a  time  dependency  of  the  form  6(t)  is 
considered.  The  appropriate  forcing  function  for  Equation  [l6] 


HYDRONAUTICS,  Incorporated 


Let  4*  be  the  solution  of  the  equation 


^  J  ^  +  Sfi  .  J  H.  N-  ,  6(^)  6(y)  S(t)  [351 


therefore,  the  vertical  velocity  component  v’  is  given  by 


v'  =  2tA  e  r—  -  b<|) 

1  ^y 


The  Fourier  Integral  representation  of  4)  is  easily  found 


to  be 


4-00  +00  -hoo 


4>(x,y,t)  =  — - — 

(2irf 


—  00  -00  —00 


1  (ax  +  j3y  +  tt)t ) 


la®  -  (j3®+  b® ) 


dadjSdo) 


Integrate  first  with  respect  to  cd.  The  Initial  condition 
of  the  fluid  is  the  state  of  hydrostatic  equilibrium.  With  this 
consideration  in  mind, the  path  of  Integration  is  displaced 
slightly  below  the  real  axis.  The  contour  is  closed  by  a  large 
semi-circle  in  the  upper  half  plane  and  the  integral  evaluated 
by  residues  to  give 


# 
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4>  = 


-.00  -00 
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The  points  of  stationary  phase  of  the  integrals  occur  at 


0, 


[38] 


As  before  we  consider  only  x  >  0.  For  x/t  >  0  only  Ii  can 
have  points  of  stationary  phase.  Given  a  pair  (x/t,  y/t)  the 
real  solutions  of  the  following  simultaneous  equations 


^  (a^  +  B®  + 


^  ^  _ No^ _ 

^  (a®  +  ]3^  +  b® 


[40] 


are  points  of  stationary  phase.  These  two  equations  may  be 
combined  to  give 


Z 

X 


i: 


+  b= 


X 


L-t 


-  1 


[41] 
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Por  a  given  y/x,  the  above  equation  passes  from  a  situation 
of  no  real  solution  to  a  situation  of  two  real  solutions  as  x/t 
tends  to  zero  from  a  very  large  value.  The  transition  where  the 
equation  has  only  one  real  solution,  therefore,  corresponds  to 
the  maximum  distance  from  the  origin  at  a  given  time  where  the 
disturbance  is  appreciable  if  we  keep  terms  of  0(l/t)  only.  Hence, 
the  wave  front  is  the  locus  of  the  transition  points  for  all  y/x. 
The  equation  of  the  wave  front  is  obtained  by  eliminating  j3  from 
Equation  [4l]  and  its  derivative  with  respect  to  The  resulting 
equation  is 


X 

t 


[42] 


where 

a®  +  V  ka*  +  3  a® 


a  *  y/x 

Since  Equation  [36]  Involves  differentiation  of  it*  with 
respect  to  y  only,  the  above  asymptotic  property  of  the  wave  front 
also  applies  to  that  of  the  internal  waves  generated  by  the  in¬ 
stantaneous  source. 

Equation  [42]  is  plotted  in  Figure  5.  It  is  seen  that  the 
wave  fronts  for  different  times  are  geometrically  similar. 
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FIGURE  I -ALLOWABLE  WAVE  NUMBERS  FOR  DIFFERENT  FREQUENCIES 
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FIGURE  2  -  DEPENDENCE  OF  THE  PHASE  VELOCITY  ON  THE  DIRECTION  OF  THE  WAVE  NUMBER 
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FIGURE  3  -  DEPENDENCE  OF  THE  MAGNITUDE  OF  THE  GROUP  VELOCITY  ON  DIRECTION  OF 

THE  WAVE  NUMBER 
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FIGURE  5-THE  WAVE  FRONT  OF  AN  INSTANTANEOUS  DISTURBANCE  IN 

AN  EXPONENTIALLY  STRATIFIED  FLUID  OF  INFINITE  EXTENT 
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